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ABSTRACT 

We present cosmological parameters derived from the angular power spectrum of the cosmic mi- 
crowave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg 2 with the At- 
acama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 
500 < I < 10000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground 
contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson 
power from radio and infrared point sources, and clustered power from infrared point sources. At 
I = 3000, about half the power at 148 GHz comes from primary CMB after masking bright radio 
sources. The power from thermal and kinetic SZ is estimated to be 63000 — 6.8 ± 2.9 /^K 2 , where 
B e = l{l+l)C e /2it. The IR Poisson power at 148 GHz is B 3000 = 7.8 ±0.7 ^K 2 {C t = 5.5 ±0.5 nK 2 ), 
and a clustered IR component is required with $3000 = 4.6 ± 0.9 /iK 2 , assuming an analytic model for 
its power spectrum shape. At 218 GHz only about 15% of the power, approximately 27 /J.K 2 , is CMB 
anisotropy at I = 3000. The remaining 85% is attributed to IR sources (approximately 50% Poisson 
and 35% clustered), with spectral index a = 3.69 ± 0.14 for flux scaling as S(v) oc v a . We estimate 
primary cosmological parameters from the less contaminated 148 GHz spectrum, marginalizing over 
SZ and source power. The ACDM cosmological model is a good fit to the data (x 2 /dof = 29/46), 
and ACDM parameters estimated from ACT+ WMAP are consistent with the 7- year WMAP limits, 
with scale invariant n s = 1 excluded at 99.7% CL (3a). A model with no CMB lensing is disfavored 
at 2.8ct. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the 
ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The 
ACT data combined with WMAP give a 6a detection of primordial helium, with Yp = 0.313 ± 0.044, 
and a 4a detection of relativistic species, assumed to be neutrinos, with N e g — 5.3 ± 1.3 (4.6 ± 0.8 
with BAO±i/o data). From the CMB alone the running of the spectral index is constrained to be 
dn s /dlnk = —0.034 ± 0.018, the limit on the tensor-to-scalar ratio is r < 0.25 (95% CL), and the 
possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension 
G>< 1.6 x 10- 7 (95% CL). 

Subject headings: cosmology: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 

Measurements of anisotropics in the cosmic microwave 
background (CMB) have played a central role in the 
development of the current concordance cosmological 
model (e.g., Smoot et al. 1992; Miller et al. 1999; de 
Bernardis et al. 2000; Hanany et al. 2000; Spergel et al. 
2003). The ACDM model describes a flat universe with 
5% normal matter, 23% dark matter, 72% dark energy, 
and power-law Gaussian primordial fluctuations consis- 
tent with simple inflationary models (see e.g., Komatsu 
et al. 2010). Its parameters have been measured to a 
few percent level accuracy, using CMB data from the 
WMAP satellite and higher resolution experiments, com- 
bined with observations of large scale structure and the 
local expansion rate (Brown et al. 2009; Riess et al. 2009; 
Reid et al. 2010; Percival et al. 2010; Larson et al. 2010; 
Komatsu et al. 2010). The model fits a range of recent 
astronomical data including Type la supernova (Hicken 
et al. 2009; Kessler et al. 2009), galaxy cluster mea- 
surements (Vikhlinin et al. 2009; Mantz et al. 2010b; 
Rozo et al. 2010) and gravitational lensing observations 
(Massey et al. 2007; Fu et al. 2008; Schrabback et al. 
2010; Suyu et al. 2010). 

The WMAP satellite has measured the CMB over the 
full sky down to 0.2° resolution. Measurements at higher 
resolution have been made by a set of complementary 
balloon and ground- based experiments (e.g., Jones et al. 
2006; Brown et al. 2009; Reichardt et al. 2009; Sievers 
et al. 2009). The Atacama Cosmology Telescope (ACT) 
now measures fluctuations on scales from ~0.4° to an 
arcminute. The signal observed in this angular range 
is composed of the damped acoustic peaks of the pri- 
mordial CMB signal (Silk 1968), subsequently lensed by 
large-scale structure, as well as point source emission 
and fluctuations due to the Sunyaev-Zel'dovich (SZ) ef- 
fect (Sunyaev & Zel'dovich 1970), in which CMB photons 
scatter off electrons in the hot intra-cluster and filamen- 
tary inter-galactic media. Limits on the SZ power spec- 
trum have been reported from the ACBAR, CBI and 
SZA experiments (Reichardt et al. 2009; Sievers et al. 
2009; Sharp et al. 2010), with a recent detection re- 
ported by the South Pole Telescope (SPT) (Lueker et al. 
2009). The first power spectrum measurement from ACT 
(Fowler et al. 2010) provided a limit on the SZ power 
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spectrum, as well as on the point source contribution. 

In this paper we present cosmological parameter con- 
straints from power spectra estimated from the ACT 
2008 observing season. We use the power spectrum to 
constrain a model for the SZ and point source contribu- 
tion in the ACT 148 GHz and 218 GHz data. We then 
combine the ACT 148 GHz data with WMAP observa- 
tions, and additional cosmological distance measures, to 
constrain the ACDM model and a set of extensions that 
have particular effects at small scales. 

This is one of a set of papers on the ACT 2008 data 
in the southern sky: Swetz et al. (2010) describes the 
ACT experiment; Dunner et al. (2010) describes the 
observing strategy and the data; Hajian et al. (2010) 
describes the calibration to WMAP; Das et al. (2010) 
presents the power spectra measured at 148 GHz and 
218 GHz, and this paper estimates parameters from the 
power spectrum results. A high-significance SZ galaxy 
cluster catalog is presented in Marriage et al. (2010a), 
with multi- wavelength observations described in Menan- 
teau et al. (2010) and their cosmological interpretation 
in Sehgal et al. (2010b). Marriage et al. (2010b) presents 
the 148 GHz point source catalog. Hincks et al. (2009) 
and Fowler et al. (2010) presented the first maps of clus- 
ters and power spectra respectively, from a preliminary 
version of these maps. Improved map-making and power 
spectrum estimates, with use of a larger area of sky, now 
allow us to place new constraints on cosmological models. 

The paper is structured as follows. In §2 we describe 
the ACT likelihood and parameter estimation method- 
ology. In §3 we present results on SZ and point source 
parameters from the small-scale power spectra. In §4 we 
present constraints on a set of cosmological models in 
combination with other cosmological data, and conclude 
in §5. 

2. METHODOLOGY 

This section describes the methods used to estimate 
parameters from the ACT power spectra. The power 
spectra, estimated from 296 deg 2 of sky observed in 2008, 
are described in Das et al. (2010), and details of the ex- 
periment, data reduction and map-making are described 
in Swetz et al. (2010) and Dunner et al. (2010). We 
will estimate two sets of parameters: 'primary' and 
'secondary'. Primary parameters describe the cosmo- 
logical model from which a theoretical primary CMB 
power spectrum can be computed. Secondary parame- 
ters describe the additional power from SZ fluctuations 
and foregrounds. We construct an ACT likelihood func- 
tion that returns the probability of the ACT data given 
some theoretical CMB power spectrum and a set of sec- 
ondary parameters. Primary and secondary parameters 
are then estimated from ACT and additional datasets us- 
ing Markov Chain Monte Carlo (MCMC) methods. The 
ACT likelihood function is described in Section 2.1, and 
the MCMC methods in 2.2. 

2.1. ACT Likelihood 

For maps of temperature fluctuations at 148 GHz and 
218 GHz, three cross-spectra are estimated in bands in 
the range 500 < I < 10000 for 148 GHz, and 1500 < 
I < 10000 for both the 218 GHz and the 148 x 218 GHz 
cross-spectrum (see Das et al. 2010). A likelihood func- 
tion is constructed to estimate parameters from these 
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spectra. The function returns the likelihood of the data, 
p(d\Cf MB , 8), given a theoretical lensed CMB spectrum, 
C™ B , and a set of secondary parameters describing the 
additional small-scale power, 9. In this analysis we con- 
sider two likelihoods: the '148+218' likelihood, which 
returns the likelihood of all three spectra given a model, 
and the '148-only' likelihood, which returns the likeli- 
hood of just the 148 GHz spectrum given a model. 

The temperature fluctuations in the ACT maps are 
expected to be the sum of fluctuations from lensed 
CMB, thermal and kinetic SZ, radio point sources, in- 
frared point sources, and thermal dust emission from the 
Galaxy (see e.g., Sehgal et al. 2010a). These vary as func- 
tions of frequency. The lensed CMB and the kinetic SZ 
are blackbody emission, and so are constant as a function 
of frequency in thermodynamic units. The thermal SZ 
emission has a known frequency dependence (Sunyaev 
& Zel'dovich 1970) and has negligible contribution at 
218 GHz. The radio point sources emit synchrotron, and 
the IR sources emit thermally, so both can be modeled 
as following power law emission in flux given a narrow 
enough frequency range. For frequency v and direction 
n the signal in the maps can be modeled in thermal units 
as 

AT(i/, n) = AT CMB (n) + AT SZ (^, n) + AT fg (^, n), (1) 

with lensed CMB fluctuations AT CMB (n). The SZ signal 
is the sum of thermal and kinetic components 



AT (v, n) 



AT tSZ (n) + AT ksz (n), (2) 



kSZ^ 



where the factor f(v) = 2 — (x/2)/ tanh(x/2), for x = 
hv/ksTcMB, converts the expected SZ emission from 
the Rayleigh- Jeans (RJ) limit to thermodynamic units, 
and ATq SZ is the expected signal at frequency v Q . At 
218 GHz there is negligible thermal SZ signal, with 
/(218) = 0. The point source and Galactic components 
are modeled as 



AT fg (^, n) : 



ATq R 



A *'"" ( .s)' 



>0 



(3) 



assuming that the IR and radio source emission in an- 
tenna temperature, ATQ R ' rad , scale with global power 
laws ad — 2 and a s — 2, respectively, where a is the index 
in flux units. The factor g(y) — (e x — l) 2 /x 2 e x converts 
from antenna to CMB thermodynamic temperature at 
frequency v. The factors are g = [1.71,3.02] for 148 GHz 
and 218 GHz. Power law behavior is also assumed for 
the Galactic dust emission AT Gal . This behavior is ex- 
pected to be a good approximation between 148 GHz and 
218 GHz, but breaks down over larger frequency ranges. 

To compute the likelihood one could first estimate the 
CMB map, by subtracting off the foreground and SZ 
components. This is commonly done in CMB analyses 
for subtracting Galactic components (e.g., Bennett et al. 
2003), and has also been used for subtracting the IR 
point sources (Hall et al. 2010; Lueker et al. 2009). How- 
ever, for the ACT frequencies and noise levels the radio 
sources cannot be neglected, so a linear combination of 



the 148 GHz and 218 GHz maps will not remove all the 
source contamination. Instead, we choose to construct a 
model for the cross power spectra between frequency Vi 
and vj, 

&> = (f;(v i )f t ( Vj )), (4) 

where is the Fourier transform of AT(fi). For ACT 
analysis we use a flat-sky approximation, described in 
Das et al. (2010). The individual components arc 
assumed to be uncorrelated, so the theoretical cross- 
spectrum Bf 1 ' 1 ' is modeled from Eq. 1 as 



B 



th,ij 



B. 



CMB 



+ S tsz,i j + S ksz, ij _ 



+ Bf Xi + B 



rad,ij _|_gGal,ij 



(5) 

where B t = 1(1 + l)C t /2-n. The first term, Bf MB , is the 
lensed primary CMB power spectrum and is the same at 
all frequencies. The thermal SZ (tSZ) power is modeled 

as 

B tsz,y 



Atsz 



/("<) /(^)gtSZ 



fiyo) fM 



J 0J ' 



(6) 



where B^f is a template power spectrum corresponding 
to the predicted tSZ emission at i/q =148 GHz for a model 
with as — 0.8, to be described in Section 2.1.2, and A t sz 
is an amplitude parameter. The kinetic SZ (kSZ) power 
is modeled as 



B 



kSZ,ij 



A RkSZ 



(7) 



where £>g^ z is a template spectrum for the predicted 
blackbody kSZ emission for a model with as = 0.8, also 
described in Section 2.1.2. The total SZ power is then 
given by 



B^ = A 



fiyj) K tsz i A K ksz 



(8) 



The infrared point sources are expected to be clustered, 
and their power is modeled as 



B 



IR,ij 



A, 



( 1 



V3000 



+ A C B* 



■clust 

£ 



u d -2 



(9) 



where Ad and A c are the values of B^ 00 at 148 GHz for 
Poisson and clustered dust terms respectively, assuming 
a normalized template spectrum £>§'" st . This will be de- 
scribed in Section 2.1.3. This model assumes the same 
spectral index for the clustered and Poisson IR power. 
The radio sources are not expected to be significantly 
clustered (see, e.g., Sharp et al. 2010; Hall et al. 2010), 
and so can be described by Poisson scale- free power, with 



R rad,ij 



A. 



3000 



g(v ) g{vo) \vo v 



(10) 



with amplitude A s normalized at vq =148 GHz and I = 
3000. 

Though we have described the pivot frequency as 
vq =148 GHz for all components in Eqs. 6-10, in practice 
we use the band-centers for SZ, radio and dusty sources 
given in Table 4 of Swetz et al. (2010). The Galactic 

emission, B e a is expected to be sub-dominant on ACT 
scales, as demonstrated in Das et al. (2010) using the 
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FDS dust map (Finkbeiner et al. 1999) as a Galactic 
dust template, so is neglected in this analysis. 

Given SZ and clustered source templates, aside from 
parameters constrained by £>f MB , this model has 7 free 
parameters: five amplitudes A t $z, A^sz, Ad, A c , A s , and 
two spectral indices, ay, a s . As we will describe in Sec 
2.2 A, wc impose priors on some of these and constrain 
others. We refer to these parameters as 'secondary' pa- 
rameters, to distinguish them from 'primary' cosmologi- 
cal parameters describing the primordial fluctuations. In 
part of the analysis we will estimate parameters from the 
148 GHz spectrum alone. In this case i = j and Vi = uq. 
The model in Eq. 5 then simplifies to 

Bf = B? MB + A tsz B*f + A kSZ B*f 

+{As + Ad) {^Oo) +AcB oT- (11) 
This can be further simplified to 

Bf = Bf MB + A SZ B^ + A v (3^5) 2 + A c^T > (I 2 ) 

where A p = A s + Ad is the total Poisson power at 
£ = 3000, and A$z = A^sz = A t $z measures the total SZ 
power, £>q^ = B^f + £>o S /- This is the same parameter- 
ization considered in Fowler et al. (2010), and has just 
three secondary parameters: A$z, A. p and A c . 

Using this model, we can compute the theoretical spec- 
tra for a given set of secondary parameters 8, and for 
a given theoretical CMB temperature power spectrum. 
The data power spectra are not measured at every mul- 
tipole, so bandpower theoretical spectra are computed 
using C^ h ' lj = wfycf 1 ' 1 ', where w % b \ is the bandpower win- 
dow function in band b for cross-spectrum ij, described 
in Das et al. (2010). 1 The data power spectra have cal- 
ibration uncertainties (to be described in Sec 2.1.1). To 
account for these uncertainties we include two calibration 
parameters y(i>i), for each map i, that scale the estimated 
data power spectra, C b , and their uncertainties, as 

<%=v{vi)v{vj)&i. (13) 
The likelihood of the calibrated data is then given by 
-2 biJgf = (C£ h - C 6 ) T S- 1 (C^ h - C b ) + lndct £, (14) 
assuming Gaussian uncertainties on the measured band- 
powers with covariance matrix S. For the 148+218 like- 
lihood the model and data vectors C b h and Cb contain 

,1 4. y-f r/~<148.148 ^-,148,218 ^218,2181 

three spectra, Cb = [C b , C b ' ,C b \. For the 
148-only likelihood, C b = C 148 ' 148 . We use the data 
between 500 < I < 10000 for the 148 x 148 GHz auto- 
spectrum, but restrict the range to 1500 < t < 10000 
for the 218 x 218 GHz and the 148 x 218 GHz spectra. 
This range is chosen since for 148 GHz at scales larger 
than t — 500 the signal cannot be accurately separated 
from atmospheric noise, and for 218 GHz the maps do 
not converge below t — 1500, described in Das et al. 
(2010). The bandpower covariance matrix S is described 
in the Appendix of Das et al. (2010), and includes cor- 
relations between the three spectra. Das et al. (2010) 
demonstrates with Monte-Carlo simulations that the co- 
variance is well modeled by a Gaussian distribution with 
negligible correlation between bands. 

1 Here we use the notation wm for the bandpower window func- 
tions; Das et al. (2010) uses B(,£. 



2.1.1. Calibration and beam uncertainty 

The ACT calibration is described in Hajian et al. 
(2010). The 148 GHz maps are calibrated using WMAP, 
resulting in a 2% map calibration error in temperature 
units, at effective I — 700. The 218 GHz maps are cal- 
ibrated using observations of Uranus, with a 7% cali- 
bration error at I = 1500. The two calibration errors 
have negligible covariance, and are treated as indepen- 
dent errors. For analyses using 148 GHz data alone we 
marginalize over the calibration uncertainty analytically 
following Ganga et al. (1996) and Bridle et al. (2002). For 
joint analyses with the 148 GHz and 218 GHz data wc 
explicitly sample the calibration parameters yiyi) with 
Gaussian priors of y(148) = 1.00 ± 0.02 and y(218) = 
1.00 ± 0.07. We check that the analytic and numeri- 
cal marginalization methods give the same results for 
148 GHz. 

The beam window functions are described in Das et al. 
(2010), and are estimated using maps of Saturn. The 
maps are made with an independent pipeline to the ini- 
tial ACT beam estimates made in Hincks ct al. (2009), 
but produce consistent results. The uncertainties on 
the beam window functions are of order 0.7% for the 
148 GHz band and 1.5% at 218 GHz. Uncertainties in the 
measured beams are incorporated using a likelihood ap- 
proximation described in Appendix A; the magnitude of 
the derived uncertainties is consistent with Hincks et al. 

(2009) and the uncertainties used in the parameter esti- 
mation in Fowler et al. (2010). 

2.1.2. SZ templates 

The thermal SZ template B^f describes the power 
from tSZ temperature fluctuations from all clusters, nor- 
malized for a universe with amplitude of matter fluctu- 
ations (Ts = 0.8. There is uncertainty in the expected 
shape and amplitude of this signal, arising due to incom- 
plete knowledge of the detailed gas physics that affects 
the integrated pressure of the clusters. Previous cosmo- 
logical studies, e.g., the ACBAR and CBI experiments, 
have used templates derived from hydrodynamical sim- 
ulations (Bond et al. 2005). The analysis for WMAP 
used the analytic Komatsu-Seljak (K-S) spectrum de- 
rived from a halo model (Komatsu & Seljak 2001). Re- 
cent studies for SPT have considered simulations and 
analytic templates from Sehgal et al. (2010a) and Shaw 
et al. (2009). 

In this analysis we consider four thermal SZ tem- 
plates, from Sehgal et al. (2010a), Trac, Bode, & Os- 
triker (2010), Battaglia et al. (2010), and Shaw et al. 

(2010) . Trac, Bode, & Ostriker (2010) constructed sev- 
eral templates by processing a dark matter simulation to 
include gas in dark matter halos and in the filamentary 
intergalactic medium. Their 'standard' model, which was 
first described in Sehgal ct al. (2010a), is referred to here 
as 'TBO-1'. It is based on the gas model in Bode, Os- 
triker, & Vikhlinin (2009), with the hot gas modeled 
with a polytropic equation of state and in hydrostatic 
equilibrium, with star formation and feedback calibrated 
against observations of local clusters. This is the tem- 
plate considered in the ACT analysis by Fowler et al. 
(2010), and has a similar amplitude and shape to the 
K-S spectrum. 

Recent state-of-the-art simulations described in 
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Battaglia et al. (2010) have been used to predict the SZ 
spectrum (referred to as 'Battaglia'). Full hydrodynam- 
ical SPH simulations were made of large scale cosmic 
structure, including radiative cooling, star formation, 
feedback from AGN and supernovae, and non-thermal 
pressure support. The predicted spectrum has 2/3 the 
power compared to the TBO-1 spectrum and is more con- 
sistent with SZ measurements from SPT (Lueker et al. 
2009). It is also compatible with predictions made that 
AGN heating would decrease the expected SZ power 
(Roychowdhury et al. 2005). 

These hydrodynamical simulations pre-date the SPT 
observations, but there have also been recent develop- 
ments in simulating and modeling the expected SZ ef- 
fect in light of the SPT results in Lueker et al. (2009), 
and motivated by recent observations of the intra-cluster 
medium (see Trac, Bode, & Ostrikcr (2010) for a dis- 
cussion). In a second model described in Trac, Bode, 
& Ostriker (2010), the nonthermal20 model referred to 
here as 'TBO-2', there is 20% non-thermal pressure sup- 
port, with increased star formation and reduced feed- 
back, which has the effect of lowering the predicted SZ 
power. The 'Shaw' model, described in Shaw et al. 
(2010), takes an analytic halo model approach, assum- 
ing hydrostatic equilibrium and a polytropic equation of 
state, with star formation, feedback from supernova and 
AGN, and energy transfer from dark matter to gas dur- 
ing mergers. It includes non-thermal pressure support. 
The spectra from each of these models are shown in Fig- 
ure 1, with all models normalized to a$ = 0.8. We study 
constraints on all four templates. Apart from the TBO-1 
template, all have similar amplitudes of S3000 ~ 5 — 6/1K 2 
at t — 3000, although the models have different amounts 
of star formation, feedback, and non-thermal support. 

The shape and amplitude of the expected kinetic SZ 
power spectrum is highly uncertain. We use the kSZ 
template described in Sehgal et al. (2010a), also shown 
in Figure 1. The corresponding template for the non- 
thermal model in Trac, Bode, & Ostriker (2010) has 
a similar amplitude. It is also consistent with predic- 
tions from second order perturbation theory (Hernandez- 
Montcagudo & Ho 2009). In this analysis the contamina- 
tion of the SZ signal by point sources is neglected, which 
is shown in Lin et al. (2009) to be a good approximation 
for radio galaxies. Lueker et al. (2009) show it is also a 
good assumption for IR sources for the current levels of 
sensitivity. 

2.1.3. Clustered source templates 

The shape and amplitude of the power spectrum of 
clustered dusty galaxies are not yet well characterized 
(Knox et al. 2001; Fernandez-Conde et al. 2008; Viero 
et al. 2009; Hall et al. 2010), although there have been 
measurements made by Viero et al. (2009) from the 
BLAST experiment. In Fowler et al. (2010) we adopted 
a simple power law model, with Bi cx I. In this analysis 
we move beyond this simple parameterization to consider 
two alternative model templates. The first, 'Src-1', is ob- 
tained from the infrared source model described in Sec- 
tion 2.5.2 of Sehgal et al. (2010a). This model assumes 
that the IR emission traces star formation in halos at 
z < 3, and that the number of IR galaxies in a given 
halo is proportional to its mass. For the spectral energy 
distribution (SED) of the galaxies, the model uses an ef- 
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Fig. 1. — Thermal SZ templates for four different models con- 
sidered in this analysis, and a single kSZ template, normalized at 
148 GHz for cosmologies with ct 8 = 0.8. The 'TBO-1' template 
is from Sehgal et al. (2010a), described further in Trac, Bode, & 
Ostrikcr (2010) together with 'TBO-2', derived from N-body simu- 
lations. The 'Battaglia' model is derived from hydrodynamic SPH 
simulations (Battaglia et al. 2010). The 'Shaw' model is based on 
an analytic halo model (Shaw et al. 2010). The 'kSZ' template is 
the kinetic SZ template in Sehgal et al. (2010a). Two clustered IR 
source templates considered ('Src-1' and 'Src-2') are described in 
Sec 2.1.3 and normalized at t = 3000. The IR source curves are 
multiplied by ten for clarity. 

fective graybody law in which the dust temperature is a 
function of the CMB temperature, but its value at z = 
is a free parameter. The dust emissivity spectral index 
and the typical IR luminosity and characteristic masses 
of the halos hosting IR galaxies are free parameters of 
the model. The model includes only the two-halo power 
spectrum given in Eq. 24 of Sehgal et al. (2010a), with 
contributions from pairs of galaxies in different halos. 
The parameters of the model have been updated from 
Sehgal et al. (2010a) to better fit the observed BLAST 
power spectra at 250-500/.tm. 2 This updated template is 
shown in Figure 1, normalized to unity at £ = 3000. The 
shape is similar to the clustered model used in the SPT 
analysis by Hall et al. (2010), peaking at i ~ 1000. 

We also consider the effect on our results of using an 
alternative template, 'Src-2', that has both one-halo and 
two-halo power, following a halo model prescription sim- 
ilar to Viero et al. (2009). Dark matter halos are popu- 
lated using galaxy source counts from the source model 
presented in Lagache, Dole, & Puget (2003), and halo 
occupation distribution parameters are tuned to fit the 
BLAST power spectra. This normalized template is also 
shown in Figure 1 . At large scales this has a similar shape 
to the Src-1 template, but at small scales tends approx- 
imately to the Bi oc I scaling adopted in Fowler et al. 
(2010), which was motivated by observations of galaxy 
clustering at small angles with typical correlation func- 
tion C{0) oc 0-°- 8 (Peebles 1980). The two templates 
differ at £ > 3000, but at these angular scales the Pois- 
son power is expected to dominate over the clustering 
term. 

2.1.4. Likelihood prescription 

2 The updated parameters are = 1.4, T = 25.5, Ah = 4X10 11 , 
M 2 = 2.5 xlO 12 , M cool = 5xl0 14 , L* = 2.3 X 10 11 , with definitions 
in Sehgal et al. (2010a). 
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To summarize the methods, an analysis with the 
'148+218' likelihood follows these steps to return the 
ACT likelihood for a given model: 

• Select primary cosmological parameters, and com- 
pute a theoretical lensed CMB power spectrum 
gCMB usm g ^ ne CAMB numerical Boltzmann code 
(Lewis et al. 2000). 

• Select values for secondary parameters = {A t $z, 
^ksz, Ad, A c , A s , ad, a s }, and compute the total 
theoretical power spectra £^ h ' lj for 148 x 148, 148 x 
218 and 218 x 218 using Eq. 5. 

• Compute the bandpower theoretical power spectra 
°6 - w M u e ■ 

• Select values for the calibration factors for 148 GHz 
and 218 GHz, and compute the likelihood using 
Eq. 14 for 500 < £ < 10000 for 148 x 148 and 
1500 < £ < 10000 for 148 x 218 and 218 x 218. 

• Add the likelihood term due to beam uncertainty, 
described in Appendix A. 

A large part of our analysis uses only the 148 GHz spec- 
trum. An analysis done with this '148-only' likelihood 
follows these steps: 

• Select primary cosmological parameters, and com- 
pute a theoretical lensed CMB power spectrum 
Bf MB using CAMB. 

• Select values for secondary parameters = {^sz, 
A p , A c } and compute the total theoretical power 
spectrum Bf 1 at 148 GHz using Eq. 5. 

• Compute the bandpower theoretical power spec- 
trum Cl h = w w Cf. 

• Compute the likelihood using Eq. 14 for 500 < £ < 
10000 for 148 GHz. 

• Add the likelihood term due to beam uncer- 
tainty, described in Appendix A, and analytically 
marginalize over the calibration uncertainty. 

2.2. Parameter estimation methods 

We use the ACT likelihood for two separate parameter 
investigations. The first uses the 148+218 likelihood to 
constrain the secondary parameters, as our initial goal is 
to characterize the small-scale behavior, and investigate 
whether this simple model sufficiently describes the ob- 
served emission. The second uses the 148-only likelihood 
to constrain primary and secondary parameters. 

2.2.1. Secondary parameters from 148 and 218 GHz 

For most of the investigation with the 148+218 likeli- 
hood we fix the primary cosmological parameters to the 
best-fit ACDM parameters estimated from WMAP, as 
our goal is to characterize the small-scale power observed 
by ACT, and check the goodness of fit of this simple 
model. To map out the probability distribution for these 
parameters we use an MCMC method. This uses the 
Metropolis algorithm to sample parameters (Metropo- 
lis et al. 1953), following the methodology described in 
Dunkley et al. (2005). 



There are seven possible secondary parameters, but we 
do not allow them all to vary freely. The radio sources de- 
tected at 148 GHz, described in Marriage et al. (2010b), 
are observed to have typical spectral index S(v) cx v~ - 5 
in flux units. By fitting a scaled source model from 
Toffolatti et al. (1998) to the detected sources, and us- 
ing it to extrapolate to fainter sources, Marriage et al. 
(2010b) predict a residual radio source power of Cg — 
2.8 ± 0.3 nK 2 . Converting units, we use these measure- 
ments to impose a Gaussian prior of A s = 4.0 ± 0.4 yuK 2 , 
and we fix a s = —0.5. We also fix A^sz = A t sz, as 
the kSZ component is subdominant at 148 GHz and the 
SZ models predict them to be the same for a given cos- 
mology. The other parameters (A t sz, Ad, A c , and a a) 
have uniform priors with positivity imposed on the am- 
plitudes. Parameter results are quoted using the means 
and 68% confidence limits of the marginalized distribu- 
tions, with 95% upper or lower limits given when the 
distribution is one-tailed. We also quote derived param- 
eters to indicate the power in different components at 
148 GHz and 218 GHz, for example the total power in SZ 
at £ = 3000, £f z 00 = (£ kSZ + B tsz ) 30 oo- 

2.2.2. Parameters from 148 GHz 

In order to explore the probability distributions for a 
set of cosmological models, we use the 148-only likeli- 
hood for parameter estimation. The focus is on using the 
1000 < £ < 3000 spectrum to improve constraints on pri- 
mary cosmological parameters. It is important that the 
SZ and foregound contribution be marginalized over, but 
we exclude the more contaminated 218 GHz data given 
the current uncertainties in the foreground model. To 
map out the distribution for cosmological parameters we 
use MCMC methods to explore the probability distribu- 
tions for various cosmological models. 

We parameterize cosmological models using 

{fl b h 2 ,n c h 2 ,n A , A^,n s ,T}. (15) 

These are the basic ACDM parameters, describing a flat 
universe with baryon density Vt[,h 2 , cold dark matter 
(CDM) density Vl c h 2 , and a cosmological constant £Ia- 
Primordial perturbations are assumed to be scalar, adia- 
batic, and Gaussian, described by a power-law with spec- 
tral tilt n s , and amplitude A^, defined at pivot scale 
fco = 0.002/Mpc. We assume 'instantaneous' reioniza- 
tion, where the universe transitions from neutral to ion- 
ized over a redshift range Az — 0.5, with optical depth r. 
Reionization likely takes place more slowly (e.g., Gnedin 
2000; Trac et al. 2008), but current CMB measurements 
are insensitive to this choice (Larson et al. 2010). We 
also consider an additional set of primary parameters 

{dn s /dlnk,r,N cS ,Y P ,Gii}, (16) 

that describe primordial perturbations with a running 
scalar spectral index dn s /d In k, a tensor contribution 
with tensor-to-scalar ratio r, a varying number of rel- 
ativistic species N e fi, varying primordial Helium frac- 
tion Yp, and cosmic strings with tension G/i, using 
the Nambu string template described in Battyc & Moss 
(2010). These parameters are added individually to the 
ACDM model in order to look for possible deviations 
from the concordance cosmology. Apart from G/i these 
parameters all take uniform priors, with positivity priors 
on r, N e g, and Yp. The tensor spectral index is fixed at 
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Fig. 2. — The angular power spectrum measured by ACT at 148 GHz and 218 GHz (Das et al. 2010), with the theoretical model for 
CMB, SZ, and point sources best-fit to the three spectra. The lensed CMB corresponds to the ACDM model with parameters derived from 
WMAP (Komatsu et al. 2010). It dominates at large scales, but falls exponentially due to Silk damping. The majority of power at i > 3000 
comes from extragalactic point sources below a Ri20 mjy flux cut after masking. The radio sources are sub-dominant, and are constrained 
by a source model fit to detected sources at 148 GHz (Marriage et al. 2010b). The infrared source emission, assumed to follow a power 
law, is dominated by Poisson power at small scale, but about 1/3 of the IR power at £ = 3000 is attributed to clustered source emission, 
assuming a template described in the text. The best-fit SZ (thermal and kinetic) contribution at 148 GHz (assuming the TBO-1 template, 
Sehgal et al. (2010a)) is 7/iK 2 at t = 3000; the subdominant kinetic SZ also contributes at 218 GHz. The data spectra and errors have 
been scaled by best-fit calibration factors of 1.02 2 , 1.02 X 1.09 and 1.09 2 for the 148 X 148, 148 X 218, and 218 X 218 spectra respectively. 



n t = —r/8, and both the index and ratio are denned as in 
e.g., Komatsu et al. (2009). The CMB power spectrum 
from cosmic strings is expected to scale as (67/j,) 2 , so we 
follow Sievers et al. (2009) and Battye & Moss (2010) 
by parameterizing the string power using q str oc (G/i) 2 . 
Limits on 67/i are then derived from g s t r - 

We generate the lensed theoretical CMB spectra using 
CAMB 3 , and for computational efficiency set the CMB 
to zero above I = 4000 where the contribution is sub- 
dominant, less than 5% of the total power. To use the 
148-only ACT likelihood there are three secondary pa- 
rameters, A$z, A p , and A c . For this part of the analysis 
we use the TBO-1 and Src-1 SZ and clustered source 
templates, checking the effect on the primary parame- 
ters of substituting alternative templates. We also im- 
pose positivity priors on these parameters. We do not 
use any information explicitly from the 218 GHz spec- 
trum in this part of the analysis, using just the 148- 
only likelihood, although results are checked using the 
148+218 likelihood. The ACT likelihood is combined 
with the seven-year WMAP data and other cosmological 
data sets. We use the MCMC code and methodology 
described in Appendix C of Dunkley et al. (2009), with 
the convergence test described in Dunkley et al. (2005). 
A subset of results are cross-checked against the publicly 

3 Version Feb 2010, with Recfast 1.5. 



available CosmoMC code. 

To place constraints on cosmological parameters we use 
the 7- year WMAP data in combination with ACT, using 
the WMAP likelihood package v4.1 described in Larson 
et al. (2010). WMAP measures the CMB over the full sky 
to 0.2° scales. All WMAP-only results shown for com- 
parison use MCMC chains from LAMBDA 4 , described 
in Larson et al. (2010). We follow the methodology 
described in Komatsu et al. (2010) to consider the ad- 
dition of distance measurements from astrophysical ob- 
servations, on the angular diameter distances measured 
from Baryon Acoustic Oscillations (BAO) at z = 0.2 and 
0.35, and on the Hubble constant. The Gaussian priors 
on the distance ratios, r s /D v (z = 0.2) = 0.1905 ±0.0061 
and r s /D v {z = 0.35) = 0.1097±0.0036, are derived from 
measurements from the Two-Degree Field Galaxy Red- 
shift Survey (2dFGRS) and the Sloan Digital Sky Survey 
Data Release 7 (SDSS DR7), using a combined analysis 
of the two data-sets by Percival et al. (2010). The param- 
eter r s is the comoving sound horizon size at the baryon 
drag epoch, and Dy{z) = [(1 + zf D 2 A (z)cz / H (z)] 1 / 3 
is the effective distance measure for angular diameter 
distance Da, and Hubble parameter H(z). The in- 
verse covariance matrix is given by Eq. 5 of Percival 
et al. (2010). The Gaussian prior on the Hubble con- 

4 http://lambda.gsfc.nasa.gov/ 
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TABLE 1 

Parameters describing SZ and extragalactic source model 
at 148 GHz and 218 GHz 



TABLE 2 

Constraints on SZ emission 



Parameter 3, 


148 + 218 GHz 


148 GHz -only 


A tS z b 
A d (mK 2 ) 
A c (pK 2 ) 
A s (pK 2 f 
a d d 


0.62 ±0.26 
7.8 ±0.7 
4.6 ±0.9 
4.1 ±0.4 

3.69 ±0.14 


< 0.77 (95% CL) 

12.0 ± 1.9 

< 7.4 (95% CL) 

4.0 ±0.4 


xVdof 


78/106 


29/46 



a The kSZ and tSZ coefficients are set equal, A^sz = A t sz- A d , A c 
and A s are the S3000 power for Poisson infrared galaxies, clustered 
infrared galaxies, and Poisson radio galaxies at 148 GHz respec- 
tively. The ACDM parameters are not varied here. 

b For the TBO-1 template. See Table 2 for other templates and 
conversion to SZ power. 

C A Gaussian prior A s = 4.0±0.4 is imposed, and index a s = —0.5 
assumed. 

d The 148 GHz-only data cannot constrain the IR point source 
index ay. 

stant, Hq — 74.2 ± 3.6 km Mpc -1 , comes from the 
magnitude-redshift relation from HST observations of 
240 low- z Type la supernovae at z < 0.1 by Riess et al. 
(2009). The error includes both statistical and system- 
atic errors. 

3. HIGH-ELL SZ AND POINT SOURCE MODEL 

In this section we determine the goodness of fit of the 
SZ and point source model described in Section 2.1 to the 
ACT 148 GHz and 218 GHz power spectra, and estimate 
its parameters. This uses the 148+218 likelihood sum- 
marized in Sec 2.1.4, initially holding the ACDM model 
fixed to the primary CMB with parameters given in Ko- 
matsu et al. (2010). The best- fit model is a good fit 
to the three ACT power spectra over the full angular 
range 500 < £ < 10000 (x 2 = 78 for 106 degrees of free- 
dom) , with constraints on parameters given in Table 1 for 
the TBO-1 SZ template and Src-1 source template. The 
spectra are shown in Figure 2, with the estimated compo- 
nents indicated at each frequency. The mean calibration 
factors, defined in Eq. 13, are 1.02 and 1.09 for 148 GHz 
and 218 GHz respectively. These are consistent with the 
expected values at the 1-1. 2a level. The best-fitting 1.09 
factor is driven by the t < 2500 part of the 148 x 218 
cross-spectrum, where the primary CMB dominates. At 
I = 3000, about half the power at 148 GHz is from the 
primary CMB (27 out of 50 ^K 2 ), with the remainder 
divided among SZ, IR Poisson and clustered power, and 
radio Poisson power (4-8 /iK 2 in each component). At 
218 GHz, only about 15% of the power comes from the 
primary CMB at t = 3000 (27 out of 170 ^K 2 ). Half 
of the power is attributed to Poisson IR sources, the re- 
maining approximately 35% to power from clustered IR 
sources. The model fits the cross-spectrum, indicating 
that a similar population of galaxies is contributing at 
both frequencies. 

3.1. Constraints on SZ power 

Using the multi-frequency spectra, power from SZ 
fluctuations is detected at more than 95% CL, with 
estimated Asz for each template (TBO-1, TBO-2, 
Battaglia, and Shaw) given in Table 2 and shown in Fig- 
ure 3, marginalized over point source parameters. The 
estimated SZ power at I — 3000 is robust to varying the 
SZ template, with total SZ power (tSZ plus kSZ) esti- 



Template 3 


A tS z h 


K SZ c 
°3000 

(pK 2 ) 


HZ. 7 

0.8 x (A^l) 


SZ.9 
0.8 x (A^ 7 ) 


TBO-1 
TBO-2 
Battaglia 
Shaw 


0.62 ±0.26 
0.96 ± 0.43 
0.85 ±0.36 
0.87 ±0.39 


6.8 ±2.9 

6.7 ±3.0 

6.8 ±2.9 
6.8 ±3.0 


0.74 ±0.05 
0.78 ±0.05 
0.77 ±0.05 
0.77 ±0.05 


0.75 ±0.04 
0.79 ± 0.04 
0.78 ± 0.04 
0.78 ± 0.04 



a Templates are from Sehgal et al. (2010a), Trac, Bode, & Ostriker 
(2010), Battaglia et al. (2010), and Shaw et al. (2010). 
b We required A^gz = A t gz, as defined in Eqs. 6-7. 
c Total tSZ and kSZ power at 148 GHz, as defined in Eq. 8. 
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Fig. 3. — One-dimensional marginalized distributions for the 
estimated thermal SZ power in the ACT power spectra. There is 
evidence at the 95% CL level for non-zero SZ power. The value 
A t gz = 1 corresponds to the predicted thermal SZ amplitude in a 
universe with ag = 0.8. The four curves correspond to the four SZ 
templates shown in Figure 1; the TBO-1 template results in a lower 
value, although all are consistent with A t gz = 1 at the 95% CL. 
The total SZ power (including kSZ) at 148 GHz and £ = 3000 for 
all the templates is consistent, with £(£ + l)Cf z /2n = 7 ± 3 fiK 2 . 

mated to be 



B 



sz 

3000 



6.8 ±2.9 fiK 2 



(17) 



The estimated template amplitude, A t sz, varies from 
0.62 ±0.26 for the TBO-1 template, to 0.96 ±0.43 for the 
TBO-2 template. Note that ^4ksz is fixed equal to A t sz 
in these cases, with amplitudes defined in Eqs. 6-7. For 
the TBO-1 template, the mean amplitude is lower than 
expected for a universe with erg = 0.8 (A t sz = 1), but 
not significantly. This is consistent with observations by 
SPT (Lueker et al. 2009), and is an improvement over 
the initial estimate of A t sz < 1-6 at 95% CL from the 
ACT power spectrum presented in Fowler et al. (2010). 
Assuming that a% is within the limits estimated from pri- 
mary CMB data, e.g. from Komatsu et al. (2010), the 
amplitude is somewhat more consistent for the TBO-2, 
Battaglia, and Shaw templates, that include more de- 
tailed gas physics, with A t sz = 1 within the 68% CL 
for these templates. In all these cases we have held the 
primary CMB parameters fixed. For a single test case 
we marginalize over the 6 primary ACDM parameters in 
addition to the secondary parameters. This marginaliza- 
tion results in an increase in the mean value of £>§qq of 
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Fig. 4. — Marginalized distributions (68% and 95% CL) for parameters describing the SZ and point source emission in the ACT power 
spectra. Left and center. The degeneracies between the total SZ power, i?f z = £(£+ l)Cf z , and the infrared point source power, £3^ R , at 
148 GHz and i = 3000 (solid unfilled contours), are broken with the addition of 218 GHz data (solid filled contours). Both the Poisson and 
clustered IR power are shown, for two different clustered source templates (solid and dashed contours). A clustered source component is 
required to fit the multi-frequency data at 5<r significance. Right: The Poisson dust power and the index a d = 3.69 ± 0.14 (power law in 
flux between 148 GHz and 218 GHz, and unconstrained from 148 GHz alone) are anti-correlated; the index indicates a dust emissivity of 
as 1.7. 



0.5 /iK 2 (a 0.2cr change), but a negligible increase in the 
uncertainty. 

The number of clusters, and therefore the expected 
SZ power, is a strong function of the amplitude of mat- 
ter fluctuations, quantified by as (Komatsu & Kitayama 
1999). In our model we scale the SZ templates by an 
overall amplitude, and would like to infer an estimate for 
as from A t sz- In Fowler et al. (2010) we assumed a sev- 
enth power scaling, with A t sz o\ (Komatsu & Seljak 
2002), giving an upper limit of erf z < 0.84 at 95% CL, for 
-Atsz < 1.6. However, the exact scaling of the shape and 
amplitude with cosmology, and in particular with as, is 
model dependent and not precisely known (Lueker et al. 
2009; Battaglia et al. 2010; Trac et al. 2010). For the 
TBO templates the combined tSZ and kSZ signal scales 
close to the 7th power, with the tSZ varying approxi- 
mately as the 8th power (Trac, Bode, & Ostriker 2010). 
To bound the possible range we compute two limits, as- 
suming the tSZ part of the template varies as either cr| 
or it|. The estimated values for af z in these cases are 
given in Table 2. 

No detections have yet been made of the kinetic SZ 
power spectrum. From SPT observations a 95% upper 
limit on B^ Q of 13/iK 2 was estimated (Hall et al. 2010). 
If we allow the kSZ amplitude to be varied independently 
of the thermal SZ amplitude, we find an upper limit from 
the ACT data on the kinetic SZ contribution of 

#30o Z o < 8 ^K 2 (95% CL). (18) 
This is consistent with predictions by Iliev et al. (2008) of 
a 2/iK 2 Ostriker- Vishniac signal at I = 3000 and a 3/iK 2 
post-reionization kSZ signal, but would exclude models 
with higher levels of kSZ from patchy reionization. 

The estimated SZ power is a small signal, less than 
10 /iK 2 , and is correlated with other parameters. We 
therefore investigate the dependence of the constraint on 
the priors imposed on other parameters. The SZ power is 
not strongly correlated with the IR point source param- 
eters when the 218 GHz data are included, as shown in 
Figure 4, and so using the 'Src-2' clustered source tem- 
plate in place of 'Src-1' has a negligible effect. There 



is some correlation with the radio source power, as this 
term also contributes predominantly at 148 GHz. Chang- 
ing the radio index to a s = has little effect, but broad- 
ening the radio prior to A s = 4±2 does reduce the signif- 
icance of the SZ detection to 0.48 ±0.27 using the TBO-1 
template; the SZ is anti-correlated with the radio power, 
and relaxing the prior allows a larger radio component 
at 148 GHz. This indicates the importance of the radio 
source characterization in Marriage et al. (2010b) for es- 
timating the SZ power at 148 GHz. A modest increase 
in the prior to A s = 4 ± 0.8 has a negligible effect. 

If we restrict the analysis to 148 GHz alone we find 
consistent results, with i3|^ 00 < 7.8 ^K 2 at 95% CL. The 
best-fit model has \ 2 = 29 for 46 degrees of freedom. 
In this case we cannot distinguish between the SZ and 
clustered source components; the joint constraint shown 
in Figure 4 shows that similar t = 3000 power limits are 
placed on both components; marginalizing over a clus- 
tered term has little effect on the estimated SZ ampli- 
tude. The mean IR Poisson term is higher in this case, 
as there are more models fitting the 148 GHz data alone 
with low SZ and clustered source power. The multi- 
frequency information then breaks this degenerecy and 
more tightly constrains the Poisson power. 

3.2. Unresolved point source emission 

The power spectrum measures fluctuations due to 
point sources below a flux cut of approximately 20 mJy. 
This is not an exact limit since the point source mask 
is constructed from sources with signal-to-noise ratio 
greater than 5 (Marriage et al. 2010b). The point source 
power observed at 148 GHz and 218 GHz has both syn- 
chrotron emission from radio galaxies, and IR emission 
from dusty galaxies. At 148 GHz the point source power, 
after removal of 5a sources, is inferred to be split in ratio 
roughly 1:2 between radio and IR galaxies. Since we im- 
pose a prior on the residual radio power from Marriage 
et al. (2010b), we do not learn new information about 
this component from the power spectrum. At 218 GHz 
the point source power is dominated by IR dust emis- 
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Fig. 5. — The power spectrum measured by ACT at 148 GHz, scaled by £ 4 , over the range dominated by primordial CMB (I < 3000). 
The spectrum is consistent with the WMAP power spectrum over the scales 500 < I < 1000, and gives a measure of the third to seventh 
acoustic peaks. The best-fit ACDM cosmological model is shown, and is a good fit to the two datasets. At I > 2000 the contribution 
from point soures and SZ becomes significant (dashed shows the total best-fit theoretical spectrum; solid is lensed CMB). Three additional 
theoretical models for the primordial CMB are shown with A r e jf=10 relativistic species, 4 He fraction Y p = 0.5, and running of the spectral 
index dn s /dlnk = —0.075. They are consistent with WMAP but are excluded at least at the 95% level by the ACT data. 



TABLE 3 

Derived constraints on unresolved IR source emission 51 





148 GHz 


218 GHz 


Poisson 


Bsooo (MK 2 ) b 


7.8 ±0.7 ±0.7 


90 ± 5 ± 10 




C<(nK 2 ) 


5.5 ±0.5 ±0.6 


63 ± 3 ± 6 




C t (JyV" 1 ) 


0.85 ±0.08 ±0.09 


14.7 ±0.7 ±1.8 


Clustered 


B3000 (mK 2 ) c 


4.6 ±0.9 ±0.6 


54 ± 12 ± 5 


Total IR 


Bsooo (MK 2 ) 


12.5 ± 1.2 


144 ± 13 



a The two errors indicate statistical uncertainty and a systematic 
error due to clustered template uncertainty. 
b Equivalent to the parameter for 148 GHz. 
c Equivalent to the parameter A c for 148 GHz. 

sion. The IR Poisson power is estimated to be Ad — 
7.8 ±0.7 fiK 2 , with derived Poisson IR power at 148 GHz 
and 218 GHz given in Table 3. A clustered component is 
required to ht the data, with A c = 4.6 ± 0.9 /iK 2 , corre- 
sponding to power at 218 GHz of 6|^ = 54 ± 12 ^K 2 . 
A model with no clustered component has a poorer fit 
to the data by A% 2 = 28, indicating a detection of clus- 
tering at the 5cr level. It is the 218 GHz power spectrum 
that provides this detection; the 148 GHz spectrum is 
consistent with no clustered component. 

In flux units, the effective index of unresolved IR emis- 
sion is 

a d = 3.69 ±0.14 (19) 
between 148 GHz and 218 GHz, where S(v) oc v a . The 



dust index and Poisson amplitude are anti-correlated, 
shown in Figure 4. This index estimate agrees with ob- 
servations by SPT, who find a = 3.9 ± 0.3 for the Pois- 
son component, and 3.8 ± 1.2 for the clustered compo- 
nent over the same frequency range (Hall et al. 2010). A 
property that can be derived from the effective index, a, 
is the dust emissivity index, /3. For galaxies at redshift 
z = the dust emission can be described by a modified 
blackbody, S(v) oc v^B v {Td)^ for dust temperature Td- 
In the Rayleigh- Jeans (RJ) limit the flux then approx- 
imates to S{u) oc vP +2 Td, with = a — 2. Using this 
relation gives a dust emissivity index measured by ACT 
of j8 = 1.7 ± 0.14, consistent with models (e.g., Draine 
2003). However, the RJ limit is not expected to be as 
good an approximation for redshifted graybodies (e.g., 
Hall et al. 2010), adding an uncertainty to (3 of up to 
~ 0.5. This should also be considered an effective in- 
dex, given the likely temperature variation within each 
galaxy. 

We test the dependence of these constraints on choices 
made in the likelihood, using the same set of tests de- 
scribed in Sec 3.1. The estimated IR source parameters 
do not depend strongly on the SZ template chosen, with 
less than O.lcr change if we use the Battaglia or TBO-1 
SZ template. If the radio source index is set to a s = 
instead of —0.5 there is a ~ 0.3<r reduction in the IR Pois- 
son source power at 148 GHz, and a 0.2cr increase in the 
spectral index. As found in Sec 3.1, if the radio source 
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Fig. 6. — One-dimensional marginalized distributions for the 6 
ACDM parameters (top two rows) derived from the ACT+ WMAP 
combination, compared to WMAP alone. The bottom row shows 

3 secondary parameters from the ACT+ WMAP data. With the 
addition of ACT data a model with n s = 1 is disfavored at the 3<r 
level. 

power uncertainty is doubled from A s — 4 ± 0.4 /iK 2 to 

4 ±0.8 /iK 2 there is only a O.lu effect. More radio source 
power can be accommodated in 148 GHz by increasing 
the width of the radio prior to 4 ± 2 /iK 2 , resulting in a 
decrease in IR Poisson power at 148 GHz of ~ la, and 
a corresponding increase in the IR index by ~ 0.8c, but 
this scenario is disfavored by the radio source counts pre- 
sented in Marriage et al. (2010b). 

Substituting the alternative halo-model 'Src-2' clus- 
tered source template reduces the estimated IR Poisson 
power by almost la. In this case the one-halo term con- 
tributes at small scales, transferring power from the Pois- 
son to the clustered component. Given our uncertainty 
in the clustered model, we adopt this difference as an 
additional systematic error on the Poisson source levels, 
shown in Table 3. In this simple model we have also as- 
sumed that the clustered and Poisson components trace 
the same populations with the same spectral index. The 
goodness of fit of the simple model supports this assump- 
tion. The detected clustering levels are compatible with 
the detections by the BLAST experiment (Viero et al. 
2009), and will be explored further in future work. 

4. COSMOLOGICAL PARAMETER CONSTRAINTS 

In this section we use the 148-only ACT likelihood to 
estimate primary cosmological parameters, in combina- 
tion with WMAP and cosmological distance priors. Fol- 
lowing the prescription in Section 2.1.4 we marginalize 
over three secondary parameters to account for SZ and 
point source contamination. We conservatively exclude 
the 218 GHz data from this part of the analysis, to avoid 
drawing conclusions that could depend on the choice of 
model for the point source power. 

4.1. The ACDM model 



The best- fit ACDM model is shown in Figure 5, using 
the combination l^Ci to highlight the acoustic peaks in 
the Silk damping regime. The estimated parameters for 
the ACT+ WMAP combination, given in Table 4 and 
shown in Figure 6, agree to within 0.5<r with the WMAP 
best-fit. The spectral index continues to lie below the 
scale invariant n s — 1, now at the 3<r level from the 
CMB alone, with n s = 0.962 ± 0.013. This supports 
the inflationary scenario for the generation of primordial 
fluctuations (Mukhanov & Chibisov 1981; Hawking 1982; 
Starobinsky 1982; Guth & Pi 1982; Bardeen, Steinhardt, 
& Turner 1983; Mukhanov, Feldman, & Brandenberger 
1992) and is possible due to the longer lever arm from 
the extended angular range probed by ACT. With the 
addition of BAO and Hq data, the significance of n s < 1 
is increased to 3.3er, with statistics given in Table 5. 

The ACDM parameters are not strongly correlated 
with the three secondary parameters (A c , A p , A$z), as 
there is limited freedom within the model to adjust the 
small-scale spectrum while still fitting the WMAP data. 
We also find consistent results if the 148+218 ACT like- 
lihood is used in place of the 148-only likelihood. 

Evidence for the gravitational lensing of the primary 
CMB signal is investigated in the companion ACT power 
spectrum paper (Das et al. 2010). A lensing param- 
eter, A^, is marginalized over that scales the lensing 
potential from Cf to A^Cf , as described in Calabrese 
et al. (2008). An unlenscd CMB spectrum would have 
Al = 0, and the standard lensing case has Al = 1. Re- 
ichardt et al. (2009) reported a detection of lensing from 
ACBAR; in Calabrese et al. (2008) this was interpreted 
as a non-zero detection of the parameter A l , with mean 
value higher than expected, Al = 3.llJ;| at 95% CL; 

Reichardt et al. (2009) estimate A L = 1.4^ at 95% 
CL from the same ACBAR data. With the ACT power 
spectrum combined with 7-year WMAP data, Das et al. 
(2010) report the measure 

A L = 1.3±g:|±i;g (68, 95% CL), (20) 
with mean value within la of the expected value. The 
goodness of fit of an unlensed CMB model has A\ 2 = 8 
worse than the best-fit lensed case, indicating a 2.8er de- 
tection of lensing. The marginalized distribution for Al 
from ACT+ WMAP, together with the standard lensed 
(Al = 1) and unlensed spectra (Al = 0), are shown 
in Das et al. (2010). The measurement adds support to 
the standard cosmological model governing the growth 
rate of matter fluctuations over cosmic time, and by ex- 
tracting information beyond the two-point function these 
measurements are expected to be improved. 

4.2. Inflationary parameters 

4.2.1. Running of the spectral index 

We constrain a possible deviation from power-law pri- 
mordial fluctuations using the running of the index, 
dn s /dhik, with curvature perturbations described by 

/ h \ Us — \n(k/ko)dn B /d In k 

Al(k) = Al(k a ) - 

(21) 

The spectral index at scale k is related to the index at 
pivot point ko by 

dn / k \ 
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TABLE 4 

ACDM AND EXTENDED MODEL PARAMETERS AND 68% CONFIDENCE INTERVALS FROM THE ACT 2008 DATA COMBINED WITH 

SEVEN-YEAR WMAP DATA. 





Parameter 3 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 








-\- dn s Jd In h 


+ r 


-1- N„a 

1 1 * Clt 


+ Y P 


+ G[i 


Primary 




2.214 ± 0.050 


2.167 ± 0.054 


2.24b ± 0.057 


2.252 ± 0.055 


2.236 ± 0.052 


2.240 ± 0.053 


A f ' T" \ A R 

AODM 


o i,2 
\l c h 


U.1127 ± U.UU54 


U.1214 ± U.UU74 


U.1U99 ± U.UUoo 


U.loz ± U.Uzo 


U.lloo ± U.UUol 


nine i n nnc tr 

U.lllo ± U.UUoo 




o , 


U. f il it U.UoU 


U.o/U it U.U40 


U. / oo ± U.UoU 


n r 79n _i_ n nqn 
U. / zU it U.UoU 


n 711 _i_ n nqi 
U. ( 11 it U.Uol 


n 7Qfi _|_ n noo. 
U. i oU it u.uzy 




7l s 


n Q^9 x n 01"? 
u.yoz zc u.uio 


i 0^9 -l- o o^q 
i.uoz it u.uoy 


n Q7A -4- o ni r 

U.yi4 It U.UIO 


n qq 1 ? -I- o 091 
u.yyo it u.uzi 


u.y i ^ it u.uio 


n Qf^Q 1 n D1 
u.yoo it u.uio 




T 


n no? _L n n 1 /i 
U.Uo f ± U.U14 


U.Uyz ± U.Ulu 


U.UoY ± U.Ulo 


U.Uoy ± U.Ulo 


n no7 _i_ n ni c; 
U.UbV ± U.Ulo 


n no7 _L n n i t: 

U.Uo r ± U.Ulo 




i r>9 A 2 

1U 


9 A 7 -1-0 11 


9 AA -1-0 11 


Z.O / I U.io 


9 AO -1 f| 19 

z.^tu m u. iz 


Z.^iO I U.ll 


9 A 1 ? -1-0 11 
Z.^iO It U. 1 1 


Extended 


dn s /d In k 
r 

N eB 

Y P 

Gfj, 




-0.034 ± 0.018 


< 0.25 


5.3 ± 1.3 


0.313 ± 0.044 


< 1.6xl0" 7 


Derived 


^8 


0.813 ±0.028 


0.841 ± 0.032 


0.803 ± 0.030 


0.906 ±0.059 


0.846 ± 0.035 


0.803 ±0.029 




Q 1)7 


0.279 ± 0.030 


0.330 ± 0.046 


0.262 ± 0.030 


0.280 ± 0.030 


0.289 ± 0.031 


0.270 ±0.029 




Ho 


69.7 ±2.5 


66.1 ± 3.0 


71.4 ±2.8 


78.9 ±5.9 


69.5 ±2.3 


70.6 ± 2.5 


Secondary 


*Smo (MK 2 ) 


< 10.2 


< 12.3 


< 10.0 


< 12.1 


< 13.0 


< 8.8 




A v ( P K 2 ) 


16.0 ±2.0 


14.9 ± 2.2 


16.0 ±2.0 


15.1 ± 2.1 


15.0 ±2.1 


16.1 ± 1.9 




A c ( M K 2 ) 


< 8.7 


< 10.4 


< 8.0 


< 11.1 


< 11.2 


< 7.4 




-2 In if 


7500.0 


7498.1 


7500.1 


7498.7 


7498.8 


7500.1 



a For one-tailed distributions, the upper 95% CL is given. For 
two-tailed distributions the 68% CL are shown. 




0.92 0.94 0.96 0.98 1.00 
n s (k Q = 0.015/Mpc) 

Fig. 7. — Two-dimensional marginalized limits (68% and 
95%) for the spectral index, n a , plotted at the pivot point 
ko = 0.015/Mpc, and the running of the index dn s /dlnk, for 
ACT+ WMAP, compared to WMAP. This model has no tensor 
fluctuations. A negative running is preferred, but the data are 
consistent with a power-law spectral index, with dn s /dlnk = 0. 

The simplest inflationary models predict that the run- 
ning of the spectral index with scale should be small (see 
e.g., Kosowsky & Turner 1995; Baumann et al. 2009), and 
the detection of a scale dependence would provide evi- 
dence for alternative models for the early universe. Cos- 
mological constraints on deviations from scale invariance 
have been considered recently by e.g., Easther & Peiris 
(2006); Kinney et al. (2006); Shafieloo & Souradeep 
(2008); Verde & Peiris (2008); Reichardt et al. (2009), 
using various parameterizations. With CMB data alone, 
the seven-year WMAP data show no evidence for sig- 
nificant running, with dn s /dhik = —0.034 ± 0.026, and 



-0.041 ±0.023 when combined with ACBAR and QUAD 
data (Komatsu et al. 2010). With the measurement of 
the power spectrum at small scales by ACT, we find 



dn s /d\nk = -0.034 ±0.018 



CL) (23) 



and dn s /d\nk =-0.024 ± 0.015 including BAO±£T . 
The estimated parameters are given in Tables 4 and 
5. Parameters are sampled using a pivot point ko — 
0.002/Mpc for the spectral index. This choice of pivot 
point results in the index being strongly anti-correlated 
with the running, with n s (0.002) = 1.032 ± 0.039. In 
Figure 7 we show the index and its running at a decor- 
related pivot point fco = 0.015/Mpc, chosen to minimize 
the correlation between the two parameters (Cortes, Lid- 
die, & Mukherjee 2007). The relation between the index 
at these two pivot points is 



n s (k 



0.015/Mpc) =n s (fc = 0.002/Mpc) 

dft 

±ln(0. 015/0.002) 



din k 



(24) 



with other cosmological parameters unchanged. The run- 
ning prefers a negative value at 1.8cr, indicating enhanced 
damping at small scales, but there is no statistically sig- 
nificant deviation from a power law spectral index. 

We choose this model to investigate the sensitivity 
of the constraints to choices made in the likelihood, as 
dn s /d\nk is more sensitive to the small-scale spectrum, 
and may be affected by the modeling of the point source 
and SZ contributions. We find less than O.lfj variation in 
primordial parameters if we substitute alternative shapes 
for the SZ and clustered source templates, or limit the 
analysis to the I < 5000 data. The beam is measured 
sufficently well over the angular range of interest that 
results are not changed if the beam uncertainty is ne- 
glected. These tests are described further in Appendix 
B, and give us confidence that the errors are not domi- 
nated by systematic effects. 
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TABLE 5 

ACDM AND EXTENDED MODEL PARAMETERS AND 68% CONFIDENCE INTERVALS FROM THE ACT 2008 DATA COMBINED WITH 

SEVEN-YEAR WMAP DATA, AND MEASUREMENTS OF Ho AND BAO 





Parameter a 




-)- dn s j d In k 


+ V 


-1- AT * 
+ "eft 


Primary 


lOOf^/i 2 


2.222 ±0.047 


2.206 ± 0.047 


2.237 ±0.048 


2.238 ±0.046 


ACDM 




0.113 ± 0.0034 


0.1148 ± 0.0039 


0.1117 ± 0.0033 


0.140 ± 0.015 




n A 


0.724 ± 0.016 


0.713 ± 0.019 


0.729 ± 0.017 


0.715 ± 0.017 




n s 


0.963 ±0.011 


1.017 ± 0.036 


0.970 ± 0.012 


0.983 ± 0.014 






U.UOU _1_ Vj.VJ-LiJ 




n ncR 4-nni5 

U. uOU _l_ \J.KJ1-\J 


n ns6 + 0014 




10 9 A|, 


2.46 ± 0.09 


2.39 ±0.10 


2.40 ±0.10 


2.44 ± 0.09 


Extended 


dn s /dink 
r 

N eS 




-0.024 ±0.015 


< 0.19 


4.56 ± 0.75 


Derived 


C8 


0.813 ±0.022 


0.820 ±0.023 


0.811 ±0.022 


0.885 ±0.039 




S~2 m 


0.276 ±0.016 


0.287 ±0.019 


0.271 ±0.017 


0.285 ±0.017 




Ho 


69.9 ± 1.4 


69.1 ± 1.5 


70.4 ± 1.5 


75.5 ±3.0 


Secondary 


B^oo (a* 2 ) 
A P (mK 2 ) 
A c (/.K 2 ) 


< 9.7 


< 11.4 


< 10.2 


< 12.1 




16.1 ± 2.0 


15.2 ± 2.0 


16.1 ± 2.0 


15.3 ±2.1 




< 8.4 


< 10.3 


< 8.4 


< 10.2 



a For one-tailed distributions, the upper 95% CL is given. For 
two-tailed distributions the 68% CL are shown. 
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Fig. 8. — Two-dimensional marginalized distribution (68% and 
95% CL) for the tensor-to-scalar ratio r, and the scalar spectral 
index n s , for ACT± WMAP data. By measuring the I > 1000 
spectrum, the longer lever arm from ACT data further breaks the 
n s — r degeneracy, giving a marginalized limit r < 0.25 (95% CL) 
from the CMB alone. The predicted values for a chaotic inflation- 
ary model with inflaton potential V(<j>) oc <f> p with 60 e-folds are 
shown for p = 3, 2, 1, 2/3; p > 3 is disfavored at > 95% CL. 

4.2.2. Gravitational waves 

The concordance ACDM model assumes purely scalar 
fluctuations. Tensor fluctuations can also be seeded at 
early times, propagating as gravitational waves. They 
contribute to the CMB temperature and polarization 
anisotropy, polarizing the CMB with both an E-mode 
and B-mode pattern (e.g., Kamionkowski et al. 1997; 
Zaldarriaga & Seljak 1997). The tensor fluctuation 
power is quantified using the tensor-to-scalar ratio r = 
^ft(^o)/A^(fco), where is the amplitude of primordial 
gravitational waves, with pivot scale k — 0.002/Mpc. 
Inflationary models predict tensor fluctuations, with am- 
plitude related to the potential of the inflaton field (see 
e.g., Baumann et al. (2009) for a recent review.) 



Direct B-mode polarization measurements from the BI- 
CEP experiment provide limits of r < 0.7 (95% CL, Chi- 
ang et al. (2010)). Temperature and E-mode fluctuations 
over a range of scales currently provide a stronger indi- 
rect constraint on r, with r < 0.36 (95% CL) from the 
WMAP data (Komatsu et al. 2010). Models with a large 
value for r have increased power at large scales, which 
can be partly compensated by increasing the spectral in- 
dex of scalar fluctuations and reducing the scalar ampli- 
tude. This L n s — r' degeneracy can be partly broken with 
lower-redshift observations that limit r < 0.24 (95% CL) 
from WMAP+BAO+iJo (Komatsu et al. 2010). It can 
also be broken by measuring temperature fluctuations 
at £ > 1000. The tightest CMB-only constraints so 
far have come from WMAP CMB data combined with 
ACBAR and QUAD small-scale CMB data (Pryke et al. 
2009), with r < 0.3 (95% CL). With ACT combined 
with WMAP we now find 

r < 0.25 (95% CL) (25) 

for the CMB temperature anisotropy power spec- 
trum alone, comparable to constraints from combined 
cosmological datasets (r < 0.19 at 95% CL for 
ACT+WMAP+BAO±ff )- The parameter estimates 
are given in Tables 4 and 5, and the dependence of 
the tensor amplitude on the spectral index is shown 
in Figure 8. For chaotic inflationary models with in- 
flaton potential V{4>) oc (ff and N e-folds of inflation, 
the predicted tensor-to-scalar ratio is r = Ap/N, with 
n s = 1 - (p + 2)/2N. The CMB data exclude p > 3 at 
more than 95% confidence for N — 60 e-folds. 

4.3. Non-standard models 

In addition to specifying the primordial perturbations, 
the concordance model assumes that there are three light 
neutrino species, that standard BBN took place with spe- 
cific predictions for primordial element abundances, and 
that there are no additional particles or fluctuations from 
components such as cosmic defects. The damping tail 
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Fig. 9. — Constraints on the effective number of relativistic species, N e g. Left: One-dimensional marginalized distribution for N B g, for 
data combinations indicated in the right panel. The standard model assumes three light neutrino species (A^g =3.04, dotted line); the mean 
value is higher, but 3.04 is within the 95% CL. Right: Two-dimensional marginalized distribution for N c h and equality redshift z cq , showing 
that N e ff can be measured separately from z cq . N c ff is bounded from above and below by combining the small-scale ACT measurements of 
the acoustic peaks with WMAP measurements. The limit is further tightened by adding BAO and Ho constraints, breaking the degeneracy 
between N c ft and the matter density by measuring the expansion rate at late times. 



measured by ACT offers a probe of possible deviations 
from this standard model. 

4.3.1. Number of relativistic species 

The CMB is sensitive to the number of relativistic 
species at decoupling. Changing the effective number 
of species affects the evolution of perturbations by alter- 
ing the expansion rate of the universe. Neutrinos also 
stream relativistically out of density fluctuations, with 
additional species suppressing the CMB peak heights and 
shifting the acoustic peak positions (Ma & Bertschinger 
1995; Bashinsky & Seljak 2004). 

The standard model of particle physics has three light 
neutrino species, consistent with measurements of the 
width of the Z boson, giving N v = 2.984 ± 0.008 (Parti- 
cle Data Book). Three neutrino species contribute about 
11% of the energy density of the universe at z ps 1100, 
with p re i = [7/8(4/ll) 4//3 A/eff]/? 7 . Cosmological datasets 
are sensitive to p re i, which can be composed of any 
light particles produced during the Big Bang that do 
not couple to electrons, ions, or photons; or any ad- 
ditional contribution to the energy density of the uni- 
verse such as gravitational waves. Three light neutrino 
species correspond to -ZV o ff=3.04. Any deviation would 
indicate either additional relativistic species, or evidence 
for non-standard interactions or non-thermal decoupling 
(Bashinsky & Seljak 2004). 

Recent constraints on the number of relativistic species 
have been explored with CMB data from WMAP com- 
bined with low redshift probes by e.g., Spergel et al. 
(2007); Ichikawa et al. (2007); Mangano et al. (2007); 
Hamann et al. (2007); Dunkley et al. (2009); Komatsu 
et al. (2010); Reid et al. (2010). With WMAP data a 
detection was made of relativistic species with N c ^> 2.7 
(95% CL), but the upper level was unconstrained. By 
combining with distance measures, Komatsu et al. (2010) 
limit the number of species to N c g— 4.34±0.88, and Reid 
et al. (2010) added optical cluster limits and LRG power 



spectrum measures to find N e g= 3.77 ± 0.67. Mantz, 
Allen, & Rapetti (2010a) include X-ray cluster gas frac- 
tion and cluster luminosity measurements from ROSAT 
and Chandra to estimate N e e= 3.4i ] 5 , improving limits 
by constraining the matter power spectrum at low red- 
shift. BBN observations limit N eS to 3.24 ±0.6 (Cyburt 
et al. 2005). 

By combining the ACT power spectrum measurement 
with WMAP, the effective number of species is estimated 
from the CMB to be 

A Gfr = 5.3 ± 1.3 (68% CL). (26) 

A universe with no neutrinos is excluded at 4c from the 
CMB alone, with the marginalized distribution shown 
in Figure 9. We can now put an upper bound on N e g 
from the CMB alone using ACT. This improved measure- 
ment comes from the third to seventh peak positions and 
heights. The right panel of Figure 9 shows the redshift 
of equality, function of the number of species. 

The relation of z eq to the number of species is given in 
Eq. 53 of Komatsu et al. (2010). With large scale mea- 
surements the observable quantity from the third peak 
height is just z eq , leading to a strong degeneracy between 
JV e ff and Q, c h 2 . With small-scale information the CMB 
data allow a measure of N e ff in addition to z cq due to 
the additional effects of anisotropic stress on the pertur- 
bations. As an example, a model with N e g= 10 that 
fits the WMAP data is shown in Figure 5. With a large 
N c s the higher peaks are damped, and slightly shifted 
to larger multipoles. The model is excluded by the ACT 
spectrum in the 1000 < I < 2500 regime. 

The central value for iV e g preferred by the 
ACT+ WMAP data is 1.7c above the concordance value, 
with increased damping over the ACDM model; im- 
proved measurements of the spectrum will help refine 
this measurement. This is not interpreted as a statisti- 
cally significant departure from the concordance value; 
the best-fit \ 2 is only 1.3 less than for A/ cff =3.04. The 



ACT 2008 Parameters 



15 



Cl 





0.7 



0.92 0.94 0.96 0.98 1.00 1.02 
n 



Fig. 10. — Constraint on the primordial helium mass fraction Yp. Left: The one-dimensional marginalized distribution for Yp derived 
from the ACT+ WMAP data compared to WMAP alone. The measurement of the Silk damping tail by ACT constrains the number of 
free electrons at recombination, giving a 6cr detection of primordial helium consistent with the BBN-predicted Yp = 0.25. Right: The 
two-dimensional marginalized distribution (68% and 95% CL) for Yp and the spectral index n s ; the degeneracy is partly broken with the 
ACT data. 



degeneracy between N c g and il c h 2 results in a higher 
mean value for cr 8 , 0.906 ± 0.059, with all estimated pa- 
rameters given in Table 4. By adding the BAO and H 
data the N c g—il c h 2 degeneracy is further broken, with 
N cS = 4.56 ±0.75 (68% CL). This central value is higher 
than from joint constraints including X-ray and opti- 
cal cluster measurements (Reid et al. 2010; Mantz et al. 
2010a); improved CMB and low redshift measurements 
will allow further constraints and consistency checks. 

4.3.2. Primordial A He abundance 

In the standard BBN model, light nuclides are synthe- 
sized in the first few minutes after the Big Bang. Mea- 
surements of the abundance of helium are therefore sen- 
sitive to the expansion rate of the universe during this 
time (Peebles 1966; Steigman et al. 1977). In standard 
BBN, the expected 4 He mass fraction, Yp, is related to 
the baryon density, flbh 2 , and the number of neutrino 
(or relativistic) species, N e g, by 

Y P = 0.2485 + 0.0016[(273.91M 2 -6) + 100(5-l)], (27) 

where S 2 = 1 + (7/43)(7V off - 3) (see e.g., Kneller & 
Steigman 2004; Steigman 2007; Simha & Steigman 2008). 
For the ACDM model, with baryon density 100fl b h 2 = 
2.214 ± 0.050 and N = 3.04 effective species, the pre- 
dicted helium fraction is Y P = 0.2486 ± 0.0006, with 
error dominated by the 0.02% uncertainty on the lin- 
ear fit in Eq. 27 (Steigman 2010). When the neu- 
trino species are allowed to vary (as in Section 4.3.1), 
the current prediction from ACT+ WMAP+B AO+H is 
Yp = 0.267 ±0.009. For comparison, the prediction from 
the baryon density derived from deuterium abundance 
measurements is Yp — 0.2482 ± 0.0007 (see Steigman 
2010, for a review). A measurement of any deviation 
from this prediction could point the way to non-standard 
models, in particular those that affect the timing of 
BBN (Steigman et al. 1977; Boesgaard & Steigman 1985; 
Jedamzik & Pospelov 2009). This includes modifications 



to the Hubble expansion rate during BBN, energy injec- 
tion due to annihilation or decay of heavy particles, parti- 
cle catalysis of BBN reactions, and time variations in fun- 
damental constants (see, e.g., Peimbert (2008); Jedamzik 
& Pospelov (2009) for discussions). 

The He abundance estimated from observations of 
metal poor extragalactic regions (see Steigman 2007; Pe- 
imbert 2008, for example) is Y P = 0.252 ± 0.004 and 
0.252 ±0.001 (Izotov et al. 2007), although a higher mea- 
surement of Y P = 0.2565±0.0010(stat)±0.0050(syst) has 
recently been made (Izotov & Thuan 2010). There are 
systematic uncertainties in the astrophysically derived 
abundances, as helium is depleted in stars. 

The CMB provides an alternative probe of the helium 
abundance when the universe was ~400,000 years old. 
Helium recombines earlier than hydrogen, at z « 1800 
rather than z « 1100, reducing the number density of 
electrons at recombination to n e = rif,(l — Yp), where n b 
is the baryon number density (Hu et al. 1995). It affects 
the CMB at small scales by modifying the recombina- 
tion process. A larger Yp decreases the electron number 
density, increasing the mean free path of Compton scat- 
tering. This leads to decreased power on small scales, 
due to enhanced Silk damping, as shown in Trotta & 
Hansen (2004) and Komatsu et al. (2010). 

For CMB analysis the primordial helium abundance 
is usually assumed to be Yp — 0.24. Constraints on a 
varying abundance from the CMB have been presented 
in Trotta & Hansen (2004); Huey et al. (2004); Ichikawa 
& Takahashi (2006); Ichikawa et al. (2008); Dunkley et al. 
(2009); Komatsu et al. (2010), with a > 3a detection of 
Y P = 0.33 ± 0.08 reported in Komatsu et al. (2010) for 
the seven-year WMAP data combined with small-scale 
CMB observations from ACBAR (Reichardt et al. 2009) 
and QUAD (Pryke et al. 2009). We now find 



Y P = 0.313 ± 0.044 



r o CL) 



(28) 



with the ACT + WMAP data combination, a significant 
detection of primordial helium from the CMB alone. The 
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Fig. 11. — Joint two-dimensional marginalized distribution (68% 
and 95% CL) for the primordial helium mass fraction Yp and the 
number of relativistic species N e ff . The two are partly degenerate, 
as increasing N c ff or Yp leads to increased damping of the power 
spectrum. The predicted standard-BBN relation between N e g and 
Yp is indicated. The concordance 7V e ff=3.04, Yp = 0.25 model 
lies on the edge of the two-dimensional 68% CL, and a model with 
N c ff=0, Yp = is excluded at high significance. 

mean value is higher than predicted from ACDM, but 
consistent at the 1.5tr level. A universe with no primor- 
dial helium is ruled out at 6cr. The distributions for Yp 
and its correlation with the spectral index are shown in 
Figure 10, with statistics in Table 4. Figure 5 shows how 
a higher helium fraction consistent with WMAP data 
(Yp = 0.5) is ruled out by ACT's determination of Silk 
damping at small scales. There is still some uncertainty 
in the exact details of recombination (e.g., Wong & Scott 
2007; Chluba et al. 2007; Switzer & Hirata 2008; Fendt 
et al. 2009; Chluba & Sunyaev 2010). A recent refine- 
ment of the numerical code for recombination used for 
this analysis (Recfast 1.5, by Seager et al. (1999), up- 
dated to match Rubino-Martin et al. (2010)), gives a 2% 
change in the spectrum at £ = 2000. This is subdomi- 
nant to the 7% percent shift from a ler change in Yp, so 
these effects are not expected to significantly affect cur- 
rent constraints, although will become more important as 
the data improve. If we consider the possible variation 
of both the primordial helium fraction and the number 
of relativistic species, the constraints on each parameter 
are weakened as there is some degeneracy between the 
two effects. The joint marginalized distribution for these 
parameters is shown in Figure 11, together with the pre- 
dicted relation between N e s and Yp assuming standard 
BBN. The concordance N eS =3.04, Y P = 0.25 model lies 
on the edge of the two-dimensional 68% CL, and a model 
with A^ff^O, Yp = 0, is excluded at high significance. 

4.3.3. Cosmic strings 

Observations of the acoustic peaks in the CMB power 
spectrum have ruled out defects from phase transitions 
as the dominant mechanism for seeding cosmic structure 
(see e.g., Vilenkin & Shellard 2000). However, certain 
inflation models predict string perturbations of similar 




1000 2000 3000 

Multipole moment I 

Fig. 12. — The power spectrum measured by ACT at 148 GHz, 
scaled by £ 4 and with best-fit secondary model subtracted, with 
the best-fit ACDM (solid) compared to a model with maximal cos- 
mic string tension Gfj, = 1.6 X 10 -7 allowed by the ACT data at 
95% CL (dashed, assuming a Nambu string template described in 
Section 4.3.3). A model with the 95% upper limit allowed without 
including ACT data, with Gfi = 2.6 X 10 — 7 , is shown for compari- 
son (dot-dashed); it overpredicts the observed power in the range 
1500 < I < 4000. 

amplitudes to the inflationary perturbations (Linde 1994; 
Dvali & Tye 1999). Using the CMB one can constrain 
the string tension, G/i, and therefore the energy scale 
at which the strings are formed. Unfortunately there is 
significant uncertainty in the predicted power spectrum 
from cosmic string-generated anisotropies, due to diffi- 
culties in modeling the string network. Most approaches 
model the network as an ensemble of string segments 
with constant average properties, with string loops pro- 
duced that decay into radiation. The equations of motion 
are solved using either the Nambu or Abelian-Higgs (AH) 
method, described in e.g., Bennett & Bouchet (1990); 
Pogosian & Vachaspati (1999); Bevis et al. (2007); Bat- 
tye & Moss (2010). 

The small-scale CMB provides a unique probe of cos- 
mic strings, with simulations and forecasts by Fraisse 
et al. (2008) and Bevis et al. (2010) predicting a power 
law behavior that could dominate over the Silk damping 
tail of the inflationary inhomogeneities, consistent with 
analytic predictions by Hindmarsh (1994). Constraints 
have been placed on the cosmic string tension from re- 
cent CMB and other cosmological data (Lo & Wright 
2005; Wyman et al. 2005, 2006; Battye et al. 2006; Fraisse 
2007; Bevis et al. 2007; Urrestilla et al. 2008; Sievers et al. 
2009). Most recently, Battye & Moss (2010) report limits 
of Gfi < 2.6 x 10~ 7 (95% CL) for Nambu strings using 
5-year WMAP data combined with large-scale structure 
and BBN data. They find a significant dependence of 
this limit on the chosen string model, with up to a factor 
of three variation. For a simple comparison, we consider 
just the Nambu cosmic string template used in Battye & 
Moss (2010), extended to scales £ > 3000 with a power 
law, Bi oc l~ Y (Fraisse et al. 2008). The template is 
held fixed for all cosmological models. Assuming this 
model, we find limits from ACT combined with WMAP 
of <7s tr < 0.025 (95% CL), which corresponds to a tension 
of 

Gpi < 1.6 x 10~ 7 (95% CL). (29) 
The spectrum corresponding to this 95% upper limit is 
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shown in Figure 12, compared to the upper limit pre- 
ACT given in Battye & Moss (2010) which overpredicts 
the power measured by ACT in the range 1500 < t < 
4000. The joint constraint on the string tension and the 
scalar spectral index may also limit the class of hybrid 
inflation models that produce cosmic strings at the end 
of inflation, and which typically predict a unity scalar 
spectral index. Battye ct al. (2006) and Bevis et al. 
(2008) demonstrated that these models provided a good 
fit to the data, but more recently Battye, Garbrecht, 
& Moss (2010) found that minimal D-term models are 
now ruled out at a 4<r level with CMB combined with 
SDSS and BBN data, and that minimal F-term models 
are increasingly disfavored. The preference from CMB 
data alone for a red spectrum with n s = 0.963 ± 0.013 
(ACT+ WMAP), marginalized over a string contribution, 
provides further evidence against these hybrid models. 

5. DISCUSSION 

The power spectra measured at 148 GHz and 218 GHz 
by ACT, using observations made in the Southern sky 
in 2008, have provided a new probe of the physics af- 
fecting microwave fluctuations at small scales. The con- 
cordance ACDM cosmological model continues to be fa- 
vored, and possible deviations from this model are more 
tightly constrained. At 148 GHz and 218 GHz, the CMB 
is dominant at scales larger than t ~ 3000 and 2000 re- 
spectively, after bright sources have been removed. At 
smaller scales, a simple model for SZ and point source 
emission is a good fit to the ACT power spectra. By us- 
ing multi-frequency information we see a preference for 
non-zero thermal Sunyaev-Zel'dovich fluctuations in the 
148 GHz power spectrum at 2a, with an amplitude con- 
sistent with observations by the South Pole Telescope. 
The mean amplitude is lower than the simplest cluster 
models predict for a universe with cr 8 = 0.8, but at less 
than 2cr significance. The level is consistent with expec- 
tations from recent models that include more complex 
gas physics; continued comparisons of observations and 
theory will allow progress to be made on cluster model- 
ing. The frequency dependence of the infrared emission 
has the expected behavior of graybody emissivity from 
dusty star-forming galaxies at redshifts 1 < z < 4, and a 
clustered infrared point source component is detected in 
the ACT data at 5a significance. 

The 1000 < I < 3000 spectrum provides a measure 
of the third to seventh acoustic peaks in the CMB, and 
the Silk damping tail from the recombination process at 
z = 1100. Using this measurement we place tighter con- 
straints on deviations from the ACDM model. Given 
the uncertainty on the expected power spectrum from in- 
frared point sources, we have limited this part of the anal- 
ysis to the power spectrum from 148 GHz. The data are 
found to be consistent with no deviations from ACDM, 
and the gravitational lensing of the temperature power 
spectrum is at the expected level, with an unlensed signal 
disfavored at 2.8<r. We have detected primordial helium 
at 6cr, and relativistic species at 4<r from the CMB alone, 
both consistent with the expected levels. The best-fit 
models prefer increased damping beyond the ACDM ex- 
pectation, but at less than the 2a level. The cosmological 
parameters considered have distinct effects on the power 
spectrum, but there is some degeneracy between N e g, 
Yp, and dn/d\nk. This means that an enhanced damp- 



ing leads to higher mean values for either N e g or Y p , or 
to more negative dn/dhik. These are seen at 1.4-1.8cr 
from the concordance value when considered individu- 
ally as extensions to ACDM, but there is no evidence 
that these deviations are statistically preferred. We do 
not find evidence for a gravitational wave component or 
a contribution from cosmic string fluctuations, indicat- 
ing the continued consistency of cosmological data with 
minimal inflationary models. 
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APPENDIX 
A. BEAM LIKELIHOOD 

Following the prescription in Hinshaw et al. (2007), the likelihood of the beam-deconvolved spectrum is written as 

C = J2(Cb - C b )(Q + Qi) 6 -J((%' -C b ,) + lndct(Q + Qi). (Al) 

where Qo contains noise and cosmic variance, and Qi contains beam error. The matrix Qi is diagonalized in the form 

Qi ss UU T , (A2) 

where U is an JVj x M matrix with M <C N b . The matrix U is computed by decomposing the unbinned beam 
covariance matrix S b , the covariance between normalized b e and bg>, as 

S b w PDP T . (A3) 

The matrix Qi is related to the beam covariance matrix H b by 

Qf=4^S b ^, (A4) 

Ol 0£ 

so the elements of matrix U are given by 

U bti = Mu[2%-Pt,iVDi\, (A5) 
where M b i are the bandpower window functions. Once decomposed, the Woodbury formula gives 

(Qo + Qi)" 1 ~(Qo + UU T )- 1 (A6) 
= Q 1 - Q X U (I + U T Q 1 U)" 1 U T Qo \ (A7) 
The likelihood is then given by C = £o + C b , where 

Co = ]T(C\ - C b )(Q ) bb XCb' -C v ) + lndet Q , (A8) 



bb' 



and 



where 



M 



C b = - - C h )M-\C v -C v ) + N M (A9) 

w 

Q^ 1 U(I + U T Q^ 1 U) _1 U T Q 1 1 (A10) 



-l 



J bb' 



and N M = lndet (I + L^Q^U). 

B. SENSITIVITY TO LIKELIHOOD ASSUMPTIONS 

A set of assumptions are made in the ACT likelihood. We choose the ACDM+running model to check their effect 
on cosmological parameters, as subtleties in the small-scale treatment can be probed more thoroughly with this 
model than with the 6-parameter ACDM. The fiducial constraints on the ACDM+running model use data between 
500 < t < 10000, including beam error and a 2% calibration error in temperature, and use the Src-1 clustered source 
template and the TBO-1 SZ template. The assumptions we test are the choice of SZ and clustered source templates, 
the dependence of results on the beam and calibration error, and the range of angular scales used. 

Substituting the alternative halo-model clustered source template, Src-2, has a negliblc effect on the primary cosmo- 
logical parameters, as shown in Figure Bl. Similarly, using the Battaglia SZ template (Battaglia et al. 2010), does not 
affect primary parameters. The distribution of the SZ amplitude is broadened in this case as the template has a lower 
amplitude (as described in Sec 3.1), and by setting an upper limit on Agz < 2 the distribution for A p is narrowed. We 
conclude from these tests that the systematic uncertainty on primary parameters from the model is small compared 
to the statistical uncertainty. 

We also test the sensitivity to the calibration and beam uncertainty. Removing the 2% calibration error in the maps 
decreases the errors on the primary parameters by up to 0.2ct. When the beam error is neglected, the constraints are 
tightened by about O.ler, but this is not a significant effect, consistent with the measurement of the beam. If instead 
the angular range is restricted to £ < 5000, the primary cosmological parameters are unaffected, but the SZ and point 
source amplitudes are more poorly determined. 
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